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ELASTOHYDRODYNAMICS OF ELLIPTICAL CONTACTS FOR MATERIALS 
OF LOW ELASTIC MODULUS* 

The work presented in the previous chapters related to ma- 
terials of high elastic modulus, like metals. In this chapter 
the analysis is extended to materials of low elastic modulus, 
like rubber. For these materials the elastic distortions are 
large, even with light loads. Another feature of the elastohy- 
drodynamics of low-elastic-modulus materials is the negligible 
effect of the relatively low pressures on the viscosity of the 
lubricating fluid. Engineering applications in which elastohy- 
drodynamic lubrication is important for low-elastic-modulus ma- 
terials include seals, human joints, tires, and a number of lu- 
bricated elastomeric-machine elements. 

The problem of line contacts, where side leakage of the 
fluid can be ignored, has been solved theoretically for low- 
elastic-modulus materials by Herrebrugh (1968), Dowson and 
Swales (1969), and Baglin.and Archard (1972). The solutions 
presented in the first two references were obtained numerically 
and are based on simultaneous solutions of the hydrodynamic and 
elasticity equations. The approximate analytical solution of 
Baglin and Archard (1972) relied on the assumption of a simpli- 
fied form for the film shape in the central region to solve the 
point-contact (ball on plane) situation. The work presented in 

*Published as Chapter 11 in Ball Bearing Lubrication by Bernard 
J. Hamrock and Duncan Dowson, John Wiley & Sons, Inc. Sept. 1981. 
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this chapter makes extensive use of the approach adopted by 
Hamrock and Dowson (1978) (1979a) and represents, to the best of 
the authors' knowledge, the first attempt at a complete 
numerical solution of the problem of isothermal elastohydro- 
dynamic lubrication of elliptical contacts for low-elastic^ 
modulus materials. No initial assumptions are made as to the 
pressure or film thickness within the contact, and lubricant 
compressibility and pressure-viscosity effects are considered. 

The basic elastohydrodynamic theory presented in Chapter 7 
has been used in which the conjunction was divided into a large 
number of equal rectangular areas and the elastic deformation at 
any point in the field was determined on the assumption that the 
pressure over each of the rectangles was constant. It was 
assumed that the solids could be considered to be semi-infinite 
bodies and hence that the normal restrictions imposed on a 
Hertzian analysis applied. The elastic deformation in the 
conjunction was represented by a double summation series of the 
pressure multiplied by an influence coefficient. 

The elastohydrodynamic analysis of low-elastic-modulus con- 
tacts presented in this chapter is intended to extend the stud- 
ies outlined in the previous chapters, but the normal safeguard 
of restricting the application to situations in which the con- 
junction dimensions are small compared with the principal radii 
of the solids should normally be observed. In addition, an 
alternative approach to the elasticity calculation would be 
necessary for situations in which thin layers of low-elastic- 

modulus material were supported on a rigid backing of 
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considerable extent, as represented by compliant surface 
bearings. 


11.1 Theoretical Formulation 


The basic theory developed in Chapter 7 is used here with 
some minor modifications. It was discovered that numerical con- 
vergence was considerably better if the dimensionless pressure 
was written as 


P = P Hz + P D 


( 11 . 1 ) 


where 

P Hz = dimensionless Hertzian pressure 

Pq = dimensionless pressure difference, representing the dif- 
ference between the hydrodynamic and Hertzian pressures 
at a given location 

By making use of equations (7.2) and (11.1) the Reynolds 
equation can be written as 

lx vfr Ixj + k 2 3Y V n 9Y/ 9X \ n 9X/ 


_1_ = 12U 
k 2 3Y y n 9Y J 


9CPH) 

3X 


( 11 . 2 ) 

Since the Hertzian pressure is known and can be obtained from 
equations (7.9) and (7.13), the solution of the Reynolds 
equation is concerned with the pressure difference Pq. In 
equation (11.2) compressibility and viscous terms are retained 
as variables even though it was mentioned earlier that the 
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effect of the relatively low pressures on viscosity and density 
is normally negligible. 

The equation for the dimensionless film thickness can be 
written as 


where 


S + 6_ + 6_ 
Hz P, 


H = H 0 + 


R-v 


(11.3) 


Hq = constant (initially guessed) 

S = geometric separation of rigid elliptical solids 


6 Hz = e1 astic deformation due to Hertzian pressure distribution 
= elastic deformation due to pressure difference 
The elastic deformation is evaluated in exactly the same way as 
outlined in Section 5.7. 

The nodal structure used in the numerical analysis that 


yielded the results presented in this chapter is shown in Figure 
11.1. This structure is different from that adopted for materi- 
als of high elastic modulus shown in Figure 7.3. Because of the 
dimensionless representation of the coordinates in Figure 11.1, 
the actual Hertzian contact ellipse becomes a circle regardless 
of the value of the ellipticity parameter k. The nodal struc- 


ture shown in Figure 11.1 was arrived at after much exploration 
in which the number of nodes in the semimajor and semiminor 


axes, as well as the distance from the center of the contact to 
the edges of the computing zone, was varied until the optimum 
arrangement was achieved. 

Apart from these points the theoretical formulation of the 
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elastohydrodynamic lubrication problem for materials of low 
elastic modulus is the same as that presented in Chapter 7, pro- 
vided that the appropriate modulus of elasticity is selected. 

11.2 Minimum- and Central-Film-Thickness Formulas 

The dimensionless grouping found to be useful in defining 
the elastohydrodynamic problem for high-elastic-modulus mater- 
ials is also valid for materials of low elastic modulus. The 
dimensionless film thickness can therefore be written as 

H = f(k,U,W,G) 

By varying each of these parameters separately and keeping the 
remaining parameters constant, as outlined in Chapter 8, the 
numerical results shown in Table 11.1 were obtained. 

The minimum film thickness shown in this table has been 
determined for 17 sets of input data from coupled solutions of 
the Reynolds and elasticity equations. From the information 
presented in Table 11.1 a least-squares-f it, minimum-film- 
thickness formula for a fully flooded, isothermal, elastohydro- 
dynamic elliptical contact for low-elastic-modulus materials can 
be written as 

Hmin = 7 *43(1 " 0.85 e~°- 31 k)u0-65w-0.21 (11.4) 

Table 11.1 gives the values for minimum film thickness ob- 
tained from the least-squares fit as defined by equation 
(11.4). The percentage difference between the minimum film 
thickness obtained from the elastohydrodynamic ell iptical- 
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contact theory H min and the minimum film thickness obtained 
from the least-squares fit equation H min i s expressed as 

E i in equation (8.13). In Table 11.1 the values of E x 
are within the range -8 to +3 percent. 

It is interesting to compare the dimensionless minimum- 
f i lm-thickness equation for materials of low elastic modulus 
equation (11.4) with the corresponding equation generated in 
Chapter 8 for materials of high elastic modulus 

H min = 3.63(1 - e -0 * 68k)u°* 68y-0. 073 q0. 49 (8.23) 

The powers of U in equations (11.4) and (8.23) are quite 
similar, but the power of W is much more significant for 
low-elastic-modulus materials. The expression showing the 
effect of the ellipticity parameter is of exponential form in 
both equations, but with quite different constants. 

A major difference between equations (11.4) and (8.23) is 
the absence of a materials parameter G in the expression for 
low-elastic-modulus materials. There are two reasons for this. 
One is the negligible effect of the relatively low pressures on 
the viscosity of the lubricating fluid, and the other is the way 
in which the role of elasticity is simply and automatically 
incorporated into the prediction of conjunction behavior through 
an increase in the size of the Hertzian contact zone 
corresponding to changes in load. As a check on the validity of 
this, case 9 of Table 11.1 was repeated with the material 
properties changed from those of nitride to those of silicone 
rubber. The results of this change are recorded as case 17 in 
Table 11.1. 
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The dimensionless minimum film thickness calculated from 
the full numerical solution to the elastohydrodynamic contact 
theory was 181.8xl0 -6 , and the dimensionless minimum film thick- 
ness predicted from equation (11.4) turned out to be 182.5x10“^. 
This clearly indicates a lack of dependence of the minimum film 
thickness for low-elastic-modulus materials on the materials 
parameter. 

There is interest in knowing the central film thickness, in 
addition to the minimum film thickness, in elastohydrodynamic con- 
tacts. The procedure used to obtain an expression for the central 
film thickness was the same as that used to obtain the minimum- 
film-thickness expression. The central-f i 1 m— th i ck ness formula for 
low-elastic-modulus materials as obtained from Hamrock and Dowson 
(1978) is 

H c = 7.3(1 - 0.72 e _0 - 28k )U°' 64 W _0 * 22 (11.5) 

A comparison of the central- and minimum-film-thickness 
equations, equations (11.5) and (11.4), reveals only slight 
differences. The ratio of minimum to central film thickness 
evident in the computed values given in Hamrock and Dowson 
(1978) ranged from 70 to 83 percent, the average being 77 
percent. 

11.3 Comparison of Different Investigators' Results 

To evaluate the dimensionless minimum-film-thickness equa- 
tion (11.4) developed by Hamrock and Dowson (1978), a comparison 
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was made between its predictions, those of the numerical 
solution obtained by Biswas and Snidle (1976), and the recent 
experimental findings of Jamison, et al. (1978). The Biswas and 
Snidle (1976) and the Jamison, et al. (1978) results are 
applicable only for k = 1. 

The Biswas and Snidle (1976) solution for the dimensionless 
film thickness can be written as 


where 


OW B - 1-96 m: 


- 0.11 

p 


( 11 . 6 ) 


«min = HminU " 0 * 5 < U - 7 > 

Mp = WU -0 ' 75 

The dimensionless groups given in equations (11.7) and (11.8) 
were first used by Moes and Bosma (1972). 

The central-film-thickness equation obtained from the 
experimental results of Jamison, et al. (1978) can be written as 


(Hclj = 2.4 Mp 0 * 075 


(11.9) 


where 


H c = H C U 


,-0.5 


( 11 . 10 ) 

Modifying equation (11.9) to provide a minimum-film-thickness 

equation by using the assumption that (H . ) _ q.78 (H ) 

min j • c 'j 

yields 


<Wj - l - 87 “j 0 ' 075 (ii.il) 
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Recall again that equations (11.6) and (11.11) are for k = 1 
only. 

From equations (11.7) and (11.8) and Table 11.1, we can 
write equation (11.4) as 

(H min ) = 8.53(1 - 0.85 e _0 - 31k )M” 0 - 21 (11.12) 

H * 

Therefore, for k = 1, equation (11.12) reduces to 

(iWl w ) H - 3- 21 ^°' 21 (11.13) 

Note that the negative exponent in equation (11.13) is 
numerically larger than the exponents in equations (11.6) and 
( 11 . 11 ). 

A comparison of the different investigators' results from 
equations (11.6), (11.11), and (11.13) is shown in Figure 11.2. 
The three equations seem to agree quite well with each other 
numerically over the range considered, but there is quite a 
discrepancy in the slopes of the lines. The Hamrock and Dowson 
(1978) prediction of film thickness is equivalent to the Biswas 
and Sn i die (1976) theoretical value at M p = 100 and to the 
Jamison, et al. (1978) experimental results at M p = 45. 
Therefore, even though the exponent on M p for the Hamrock 
and Dowson (1978) results is larger than those obtained by 
Biswas and Snidle (1976) and by Jamison, et al. (1978), the 
agreement is quite good for k = 1. 

Also shown in Figure 11.2 is the rigid isoviscous solution 
obtained from Kapitza (1955). The enormous difference in slopes 
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between the Kapitza (rigid) line and the three sets of results 
presented by equations ( 11.6), (11.11), and (11.13) 
demonstrates the tremendous potential of elastohydrodynamic 
action for the generation and preservation of satisfactory 
lubricating films under severe operating conditions. The 
discrepancies between the predictions of the three relationships 
represented by equations (11.6), (11.11), and (11.13) are seen 
in better perspective if they are compared with the predictions 
of the Kapitza (rigid) theory at values of M p of 100 and 
1000. 

The variation of the ratio H min /H mi is shown in Figure 
11.3, where r is the minimum film thickness for rectangular 

contacts, with the el 1 ipticity parameter k for both high- and low- 
elastic-modulus materials. If it is assumed that the minimium film 
thickness obtained from the elastohydrodynamic analysis of 
elliptical contacts can only be obtained to an accuracy of 
3 percent, we find that the ratio H min /H mi n> r approaches the 
limiting value of unity at k = 5 for high-elastic-modulus ma- 
terials. For low-elastic-modulus materials the ratio approaches 
the limiting value of unity more slowly, but it is reasonable to 
state that the rectangular-contact solution will give a very 
good prediction of the minimum film thickness for conjunctions 
in which k exceeds about 11. 


10 


11.4 Contour Plots of Results 

Contour plots of the dimensionless pressure are shown in 
Figure 11.4 for two extreme values of the dimensionless speed 
parameter U of 0.05139x10“^ and 0.5139x10“^. As mentioned 
in Chapters 8 and 9, the + symbol indicates the center of the 
Hertzian contact in each case. Because of the dimensionless 
representation of the X and Y coordinates the actual 
Hertzian contact ellipse becomes a circle regardless of the val- 
ue of k. The Hertzian contact circle is shown in each figure 
by asterisks. On each figure the contour labels and each 
corresponding value of the dimensionless pressure are given. 

The inlet region is to the left and the exit region to the right 
in each figure. 

The pressure contours shown in Figure 11.4 are nearly cir- 
cular or Hertzian. In Figure 11.4(b), the high-speed case, the 
pressure at any point in the inlet is higher than in the low- 
speed case shown in Figure 11.4(a). Inside the Hertzian contact 
region the contour values of the dimensionless pressure for the 
low- speed case are higher than those for the high-speed case. 

The pressure spikes found when dealing with materials of high 
elastic modulus (Chapter 8) are not evident in these solutions 
for low-elastic-modulus materials. The absence of a pressure 
spike for low-elastic-modulus materials has been noted before 
for nominal line contacts and is due to the pressures generated 
for a given load in a contact between low-elastic-modulus 
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materials being considerably lower than those generated in a 
contact between high-elastic-modulus materials. 

Contour diagrams of the dimensionless film thickness are 
shown in Figure 11.5 for these same values of U = 0.05139xl0 -7 
and 0.5139xlCT 7 . Figure 11.5(a) shows three regions of minimum 
film thickness: one close to the rear edge of the Hertzian 

ellipse, and two off to the side. At the higher speed 
(U = 0.5139xl0 -7 ) the minimum-film-thickness region lies on 
the midplane of the contact in the direction of rolling, between 
the center of the contact and the trailing edge of the Hertzian 
ell ipse. 

The variation of pressure and film thickness in the direc- 
tion of rolling along a line near the midplane of the conjunc- 
tion is shown in Figure 11.6 for the same two values of U con- 
sidered in Figures 11.4 and 11.5. For all the solutions ob- 
tained at various speeds, the values of the dimensionless load, 
materials, and ellipticity parameters were held fixed. Figure 
11.6 shows that the pressure at any point in the inlet region 
increases as the speed increases. The dominant effect of the 
dimensionless speed parameter on the minimum film thickness in 
elastohydrodynamic contacts for low-elastic-modulus materials 
evident in equation (11.4) is reflected in Figure 11.6. Similar 
results were found for high-elastic-modulus materials, as noted 
in Chapter 8. 

Contour plots of dimensionless pressure for the two extreme 
values of the dimensionless load parameter W that were inves- 
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tigated, 0.2202xlCT 3 and 2.202xl0 -3 , are shown in Figure 
11.7. Again the pressure contours are nearly circular, or 
Hertzian. 

Contour plots of dimensionless film thickness for the same 
two values of W are shown in Figure 11.8. In Figure 11.8(a), 
for the low-load case (W = 0.2203xl0 -3 ) , the minimum film thick- 
ness occurs directly behind the center of the contact. Likewise 
in Figure 11.8(b), for the high-load case (W = 2.202xlCT 3 ), the 
minimum film thickness also occurs directly behind the center of 
the contact but closer to the Hertzian circle. The two contours 
marked C in Figure 11.8(b) indicate a slight increase in film 
thickness before the minimum-film-thickness region is reached. 

The variation of pressure and film thickness in the rolling 
direction along a line close to the midplane of the conjunction is 
shown in Figure 11.9 for the same two values of the dimension- 
less load parameter. Once again the essential features of the 
minimum-film-thickness equation are reflected in this figure since 
a change in the dimensionless load parameter of one order of mag- 
nitude produces a considerable change in the pressure profile but 
not such a significant change in the film thickness. The small 
effect of load on minimum film thickness is a reasonable and most 
important feature of elastohydrodynamic lubrication. 

11.5 Effect of Lubricant Starvation 

By using the theory and numerical procedures outlined earlier 
in this chapter, we can investigate the influence of lubricant 
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starvation on minimum film thickness in elliptical elastohydrody- 
namic conjunctions formed by low-elastic-modulus materials. Lu- 
bricant starvation is studied by simply moving the inlet boundary 
closer to the center of the conjunction, as described in Chapter 9. 

Table 11.2 shows how the dimensionless inlet distance affects 
the dimensionless film thickness for three groups of dimensionless 
load and speed parameters. For all the results presented in this 
section the dimensionless materials parameter 6 was fixed at 
0.4276, and the ellipticity parameter k was fixed at 6. The 
results shown in Table 11.2 clearly indicate the adverse effect 
of lubricant starvation in the sense that, as the dimensionless 
inlet distance m decreases, the dimensionless minimum film 
thickness H m - n also decreases. 

Table 11.3 shows how the three groups of dimensionless speed 
and load parameters affect the limiting location of the dimension- 
less critical inlet boundary distance m* at which starvation 
starts to influence film thickness. Also given in this table are 
corresponding values of the dimensionless minimum film thickness 
for the fully flooded condition, as obtained by interpolating the 
numerical values. By making use of Table 11.2 and following the 
procedure outlined in Chapter 9, we can write the critical dimen- 
sionless inlet boundary distance at which starvation becomes im- 
portant for low-elastic-modulus materials as 


0.16 


m * - 1 + l - 0? [(t ) 2 


(11.14) 


where . is obtained from equation (11.4). 
min 

Table 11.4 shows how m* affects the ratio of minimum film 
thickness in the starved and fully flooded conditions ^ m i n> s^min # 
The dimensionless minimum film thickness for a starved condition for 
low-elastic-modulus materials can thus be written as 



(11.15) 


Therefore, whenever m < m*, where m* is defined by 
equation (11.14), a lubricant starvation condition exists. When 
this is true, the dimensionless minimum film thickness is 
expressed by equation (11.15). If m > m*. a fully flooded 
condition exists and equation (11.4) can be used to predict the 


minimum film thickness. 

In Figure 11.10 contour plots of the dimensionless pressure 
(P = p / E ' ) are shown for the group 3 conditions recorded in 
Table 11.2 and for dimensiosnless inlet distances of 1.967, 
1.333, and 1.033. Note that the contour levels and intervals 
are identical in all parts of Figure 11.10. In Figure 11.10(a), 
with m= 1.967, an essentially fully flooded condition exists. 
The contours are almost circular and extend further into the 
inlet region than into the exit region. In Figure 11.10(b), 
with m = 1.333, starvation is influencing the distribution of 
pressure and the inlet contours are slightly less circular than 
those shown in Figure 11.10(a). By the time m falls to 1.033 
(Figure 11.10(c)) the conjunction is quite severely starved and 
the inlet contours are even less circular. 
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In Figure 11.11 contour plots of the dimensionless film 
thickness (H = h/R x ) are shown, also for the group 3 condi- 
tions recorded in Table 11.2 and for dimensionless inlet dis- 
tances of 1.967, 1.333, and 1.033. These film thickness con- 
tours correspond to the pressure results shown in Figure 11.10. 
The central portion of the film thickness contours becomes more 
parallel as starvation increases and the minimum-f i lm-thickness 
area moves to the exit region. The values of the film thickness 
contours for the most starved condition (Figure 11.11(c)) are 
much lower than those for the fully flooded condition (Figure 
11.11(a)). 

Figure 11.12 more clearly describes these film thickness 
results. It shows the variation of the dimensionless film 
thickness in the rolling direction for four values of the di- 
mensionless inlet distance. The value of Y was held fixed 
near the axial center of the contact, and the dimensionless pa- 
rameters U and W were held constant as shown in group 3 of 
Table 11.2 for these calculations. This figure clearly shows 
that the central region becomes flatter as starvation becomes 
more severe. Also, in going from a fully flooded condition to a 
starved condition the film thickness decreases substantially and 
the location of the minimum film thickness moves closer to the 
exit region. 

The variation of the dimensionless film thickness perpen- 
dicular to the rolling direction is shown in Figure 11.13 for 
four values of the dimensionless inlet distance. The value of 
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X was held constant near the axial center of the contact, and 
the dimensionless parameters U and W were held constant at 
the values recorded for the conditions corresponding to group 3 
of Table 11.2. The results shown in this figure clearly indi- 
cate that the central region of the conjunction is quite flat in 
the fully flooded situation but that the curvature of the pro- 
file increases as the severity of starvation increases. 

11.6 Closure 

By modifying the procedures outlined in Chapter 7 we have 
investigated the influence of the ellipticity parameter k and 
the dimensionless speed U, load W, and materials G parame- 
ters on minimum film thickness for materials of low elastic mod- 
ulus. The ellipticity parameter was varied from 1 (a ball-on- 
plane configuration) to 12 (a configuration approaching a line 
contact). The dimensionless speed and load parameters were each 
varied by one order of magnitude. Seventeen cases were used to 
generate the minimum- and central-film-thickness relations: 

Hjnin = 7.43(1 - 0.85 e -0 ' 31k )U 0 ' 65 W -0 * 21 
H c = 7.32(1 - 0.72 e-0-28kju0.64y-0.22 

The influence of lubricant starvation on minimum film 
thickness in starved elliptical elastohydrodynamic conjunctions 
has also been investigated for materials of low elastic modu- 
lus. Lubricant starvation was studied simply by moving the 
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inlet boundary closer to the center of the conjunction in the 
numerical solutions. The results show that the location of the 
dimensionless critical inlet boundary distance m* denoting the 
location between the fully flooded and starved conditions can be 
expressed simply as 



That is, for a dimensionless inlet distance m less than m*, 
starvation occurs; and for m 2 m*, a fully flooded condition 
exists. Furthermore it has been possible to express the minimum 
film thickness for a starved condition as 


Bain, 1 




0.22 


Contour plots of pressure and film thickness in and around 
the contact have been presented for both fully flooded and 
starved lubrication conditions. It is evident from these fig- 
ures that the inlet pressure contours become less circular and 
closer to the edge of the Hertzian contact zone and that the 
film thickness decreases substantially as the severity of star- 
vation increases. 

The results presented in this chapter have revealed the 
essential features of both fully flooded and starved, ellipti- 
cal, elastohydrodynamic conjunctions for materials of low elas- 
tic modulus. 


18 



SYMBOLS 


A 

A*, B* t C*, 
D*, L*, M* , 


a 


a 


B 

b 

b 


C 


C^» • • • ,c 8 

c 

c 


D 

D 

D 

De 

d 

d 


d ; 


d i 


constant used in equation (3.113) 
relaxation coefficients 

2 

drag area of ball, m 

semimajor axis of contact ellipse, m 

a/2m 

total conformity of bearing 
semiminor axis of contact ellipse, m 
b/2m 

dynamic load capacity, N 
drag coefficient 
constants 

19,609 N/cm 2 (28,440 lbf/in 2 ) 

number of equal divisions of semimajor axis 

distance between race curvature centers, m 

material factor 

defined by equation (5.63) 

Deborah number 
ball diameter, m 

number of divisions in semiminor axis 
overall diameter of bearing (Figure 2.13), m 
bore diameter, m 
pitch diameter, m 

pitch diameter after dynamic effects have acted on ball, m 
inner-race diameter, m 
outer-race diameter, m 
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E 


E‘ 


E 

E 

<? 

I 


e 

F 

F* 

#v 

F 

F 



G 

G* 

G 


g 


modulus of elasticity, N/m 


2 


effective elastic modulus, 




? 2 

internal energy, m /s 
processing factor 

ttfmin * H n,in)/ H min3 x WO 

2 1 / <- 

elliptic integral of second kind with modulus (1 - 1/k ) 

approximate elliptic integral of second kind 

dispersion exponent 

normal applied load, N 

normal applied load per unit length, N/m 

lubrication factor 

integrated normal applied load, N 

centrifugal force, N 

maximum normal applied load (at = 0), N 


applied radial load, N 


applied thrust load, N 

normal applied load at angle f, N 

o 1/2 

elliptic integral of first kind with modulus (1 - 1/k^) 

approximate elliptic integral of first kind 

race conformity ratio 

rms surface finish of ball, m 

rms surface finish of race, m 

dimensionless materials parameter, aE 

fluid shear modulus, N/m^ 

hardness factor 

p 

gravitational constant, m/s 
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H 

A 

H 

H 


c,s 


H min 


U 

min,r 


LI 

min,s 


H min 


H min 


h 

h 


h. 

1 


8/3 2 

dimensionless elasticity parameter, W / U 

3 2 

dimensionless viscosity parameter, GW / U 
dimensionless film thickness, h/R x 
dimensionless film thickness, H(W/U)^ = F h/u hqR x 
dimensionless central film thickness, h c /R x 
dimensionless central film thickness for starved 
lubrication condition 
frictional heat, N m/s 

dimensionless minimum film thickness obtained from EHL 
elliptical-contact theory 

dimensionless minimum film thickness for a rectangular 
contact 

dimensionless minimum film thickness for starved 
lubrication condition 

dimensionless central film thickness obtained from 
least-squares fit of data 

dimensionless minimum film thickness obtained from 
least-squares fit of data 

dimensionless central-film-thickness - speed parameter, 

H C U “ 0 ' 5 

dimensionless minimum-film-thickness - speed parameter, 
H min u ~°* 5 

new estimate of constant in film thickness equation 
film thickness, m 
central film thickness, m 
inlet film thickness, m 
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h 


m 



J* 

J 



L 


50 


l 

I 

M 


film thickness at point of maximum pressure, where 
dp/dx = 0, m 

minimum film thickness, m 
constant, m 

diametral interference, m 

2 

ball mass moment of inertia, m N s 
integral defined by equation (3.76) 
integral defined by equation (3.75) 
function of k defined by equation (3.8) 
mechanical equivalent of heat 

2 

polar moment of inertia, m N s 
load-deflection constant 
el 1 ipticity parameter, a/b 
approximate ellipticity parameter 
thermal conductivity, N/s °C 
lubricant thermal conductivity, N/s °C 
fatigue life 
adjusted fatigue life 

reduced hydrodynamic lift, from equation (6.21) 
lengths defined in Figure 3.11, m 

fatigue life where 90 percent of bearing population will 
endure 

fatigue life where 50 percent of bearing population will 
endure 

bearing length, m 

constant used to determine width of side-leakage region 
moment, Nm 


22 


M 


M 


g 

p 


M 


s 


m 


m* 


m 

m 


N 

n 

n* 

"n 

P 

P 


D 


Hz 


max 


iv,as 


m 


gyroscopic moment, Nm 

-0 75 

dimensionless load-speed parameter, WU 

torque required to produce spin, N m 
2 

mass of nail, N s / m 

dimensionless inlet distance at boundary between fully 
flooded and starved conditions 
dimensionless inlet distance (Figures 7.1 and 9.1) 
number of divisions of semimajor. or semiminor axis 
dimensionless inlet distance boundary as obtained from 
Wedeven, et al. (1971) 
rotational speed, rpm 
number of balls 
refractive index 

constant used to determine lengtn of outlet region 

dimensionless pressure 

dimensionless pressure difference 

diametral clearance, m 

free endplay, m 

2 

dimensionless Hertzian pressure, N/m 

2 

pressure, N/m 

2 

maximum pressure within contact, 3F/2nab, N/m 

2 

isoviscous asymptotic pressure, N/m 
solution to homogeneous Reynolds equation 
thermal loading parameter 

2 

dimensionless mass flow rate per unit width, qnQ/pQE'R 
reduced pressure parameter 

2 

volume flow rate per unit width in x direction, m / s 
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2 

volume flow rate per unit width in y direction, m / s 
curvature sum, m 

arithmetical mean deviation defined in equation (4.1), 
operational hardness of bearing material 
effective radius in x direction, m 
effective radius in y direction, m 
race curvature radius, m 

radii of curvature, m 

cylindrical polar coordinates 

spherical polar coordinates 

defined in Figure 5.4 

geometric separation, m 

geometric separation for line contact, m 

empirical constant 

shoulder height, m 

T 0^ p max 

tangential (traction) force, N 
temperature, °C 
ball surface temperature, °C 
average lubricant temperature, °C 
ball surface temperature rise, °C 

<V p max > k=1 
viscous drag force, N 

time, s 

auxiliary parameter 

velocity of ball-race contact, m/s 



u 

u 

71 

AU 

V 

W 

w 

X 

Y 


x t> V t 

x,. Y a , z ; 


z 

Z 1 

~ ^ 

x, x, x, x 1 

y, y, y, 

z, z, 7, 7 1 

a 


p f 

e' 

r . 

Y 

Y 


velocity of ball center, m/s 

dimensionless speed parameter, tiqU/E'R x 

surface velocity in direction of motion, (u^ + u b )/2, m/s 

number of stress cycles per revolution 

sliding velocity, u & - u b , m/s 

surface velocity in transverse direction, m/s 

2 

dimensionless load parameter, F/E'R 

surface velocity in direction of film, m/s 

dimensionless coordinate, x/R^ 

dimensionless coordinate, y/R x 

dimensionless grouping from equation (6.14) 

external forces, N 

constant defined by equation (3.48) 

viscosity pressure index, a dimensionless constant 

coordinate system 

2 

pressure-viscosity coefficient of lubrication, m / N 

radius ratio, Ry/R x 

contact angle, rad 

free or initial contact angle, rad 

iterated value of contact angle, rad 

curvature difference 

viscous dissipation, N/m s 

total strain rate, s - ^ 

elastic strain rate, s - ^ 

viscous strain rate, s" * 
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£ 


G 1 

e 2 

K, 


? a 

n 

n 

n 0 

^oo 

Q 

A 


X 


flow angle, deg 

total elastic deformation, m 

lubricant viscosity temperature coefficient, °C * 
elastic deformation due to pressure difference, m 
radial displacement, m 
axial displacement, m 
displacement at some location x, m 
approximate elastic deformation, m 
elastic deformation of rectangular area, m 
coefficient of determination 
strain in axial direction 
strain in transverse direction 
angle between ball rotational axis and bearing 
centerline (Figure 3.10) 
probability of survival 

2 

absolute viscosity at gauge pressure, N s/m 
dimensionless viscosity, ti/hq 

p 

viscosity at atmospheric pressure, N s/m 
6.31xl0~ 5 N s/m 2 (0.0631 cP) 
angle used to define shoulder height 
film parameter (ratio of film thickness to composite 
surface roughness) 

equals 1 for outer-race control ano 0 for inner-race 
control 

second coefficient of viscosity 

Archard-Cowking side-leakage factor, (1 + 2/3 <* )”* 

relaxation factor 
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p 0 
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°1 

T 

T o 
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t L 
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*1 

♦ 

°0 

0) 

“B 

u b 


coefficient of sliding friction 
7/n 

Poisson's ratio 

divergence of velocity vector, ( au/ax) + (av/ay) + (aw/az), s 1 

2 , 4 

lubricant density, N s /m 
dimensionless density, p/p q 

2 4 

density at atmospheric pressure, N s /m 
2 

normal stress, N/m 

2 

stress in axial direction, N/m 
shear stress, N/m 

2 

maximum subsurface shear stress, N/m 
2 

shear stress, N/m 

2 

equivalent stress, N/m 
limiting shear stress, N/m 2 

ratio of depth of maximum shear stress to semiminor axis of 
contact ellipse 

pH 3/2 

<*>k.l 

auxiliary angle 
tnermal reduction factor 
angular location 
limiting value of 'P 

absolute angular velocity of inner race, r ad/s 
absolute angular velocity of outer race, rad/s 
angular velocity, rad/s 

angular velocity of ball— race contact, rad/s 
angular velocity of ball about its own center, rad/s 
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Subscripts: 

a 

b 

c 

be 

IE 

IR 

i 

K 

min 

n 


o 

PVE 

PVR 

r 

s 

x,y,z 

Superscript: 
( — ) 


angular velocity of ball around shaft center, rad/s 
ball spin rotational velocity, rad/s 

solid a 
solid b 
central 
ball center 

isoviscous-elastic regime 

isovi scous-rigid regime 

inner race 

Kapitza 

minimum 

iteration 

outer race 

piezovi scous-elastic regime 
piezovi scous-rigid regime 
for rectangular area 
for starved conditions 
coordinate system 

approximate 
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Figure 11. 1 - Nodal structure used for fully flooded numerical calculations. 
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Figure 11.2. - Comparison of different investigators' results. 
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Figure 11. 3. - Effect of ellipticity parameter on 
ratio of dimensionless minimum film thick- 
ness for rectangular contact, for EHL high- 
and low-elastic-modulus analyses. 
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Figure 11.6. - Variation of dimensionless pressure and film thickness on X axis for 
dimensionless speed parameters U of 0. 05139x1 0" 7 and 0.5139x10*'. The value of Y is 
held fixed near axial center of contact 
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Figure 11. 13. - Variation of dimensionless film thickness on 
Y axis for four dimensionless inlet distances and for group 3 
of Table 11. 2. The value of X is held fixed near axial center 
of contact. 
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